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O ; Abstract 

^ We undertake the study of the charged Higgs bosons predicted by the model with 

gauge symmetry SU (3) c ® SU (3)l ® U (l)x- By considering Yukawa mixing couplings 
between small (~ GeV) and large (~ TeV) scales, we show that the hypercharge-one 
and hypercharge-two H 2 Higgs bosons predicted by the model, can be simulta- 
neously produced in pp collisions at different production rates. At low energy, the 
bosons exhibit the same properties as the charged Higgs bosons from a two Higgs 
doublet model (2HDM), while H 2 are additional like-charged Higgs bosons from the 
underlying 3-3-1 model. Thus, the identification of multiple like-charged Higgs boson 
resonances may test the compatibility of theoretical models with experimental data. 
We study Hf 2 pair and associated tbH^ 2 productions at CERN LHC collider. In 
particular, we obtain that pair production can be as large as the single production 
in gluon-gluon collisions due to the interchange of a heavy neutral Z' gauge boson 
predicted by the model. By considering decays to leptons Hf 2 — > tv t , we obtain sce- 
narios where small peaks of i/^-boson events in transverse mass distributions can be 
identified over the background. 



1 Introduction 



The mechanism that breaks the electroweak symmetry, which provides masses to all elemen- 
tary particles is a subject of high priority in particle physics. In the Standard Model (SM) 
[I], the inclusion of one SU(2)l scalar doublet leads to the symmetry breaking and gener- 
ates the masses of the fundamental particles at the scale of the vacuum expectation value 
v = 246 GeV. One of the consequences of this mechanism is the prediction of one massive 
neutral scalar boson: the SM Higgs Boson (SMHB). Searches for Higgs bosons have been 
carried out by different experiments at CERN-LEP [2], Fermilab-Tevatron and CERN-LHC 
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colliders. Recently, the ATLAS and CMS collaborations at LHC obtained evidence of an 
scalar resonance around the 126 GeV [3j, [1]. However, further analysis will be necessary to 
test the compatibility of this signal with the SMHB. Other searches have excluded the pro- 
duction of neutral Higgs bosons in other regions of the Higgs mass. At Fermilab-Tevatron, 
the region 156 GeV < m H < 177 GeV [5J have been excluded. At CERN-LHC, the ATLAS 
collaboration have excluded the regions, 131 — 238 GeV, and 251 — 466 GeV j6], while the 
CMS collaboration have excluded the mass range 127 — 600 GeV at 95 % CL [7j. 

Even if the existence of the SMHB is confirmed by the future experimental analysis, it 
seems "natural" that the SM must be extended to describe physics up to the Planck scale 
(2.4 x 10 18 GeV). On one hand, the mass of the SMHB is expected to be at the electroweak 
scale, as required from unitarity constraints [8j. However, on the other hand, if the limit 
of validity of the SM is near the Planck scale, there will arise huge quantum corrections to 
the squared Higgs mass, and "unnatural" cancellations must occur to give a Higgs boson 
of the order of the electroweak scale [9j. Thus, it is to expect some sort of theory beyond 
the SM that fill the enormous range between the electroweak scale and the Planck scale. In 
particular, many theoretical extensions of the SM, charged Higgs bosons, H , are predicted. 
The detection of a charged Higgs boson could reveal many features about the underlying 
model beyond the SM. Searches at LEP in a general Two Higgs Doublet Model (THDM) 
[TU] have set a limit M H ± > 79.3 GeV [UJ. For light charged Higgs bosons (M H ± < m top ), 
in the framework of the Minimal Standard Supersymmetric Model (MSSM) [12], the main 
production channel at the CERN-LHC is through top-quark decays t — > bH + , while for 
heavy Higgs (M H ± > m top ), the associated tb production pp — > tbH ± + X is the dominant 
mode [T3] . 

Other interesting alternative to extend the SM are the models with gauge symmetry 
SU(3) C eg) SU(3)l <S> U(1)x, also called 3-3-1 models, which introduce a family nonuniversal 
£7(1) symmetry JHJ [151 EH E]- These models have a number of phenomenological advan- 
tages. First of all, from the cancellation of chiral anomalies [18] and asymptotic freedom in 
QCD, the 3-3-1 models can explain why there are three fermion families. Secondly, since the 
third family is treated under a different representation, the large mass difference between 
the heaviest quark family and the two lighter ones may be understood [19]. Finally, these 
models contain a natural Peccei-Quinn symmetry, necessary to solve the strong-CP problem 
[2D]. 

The 3-3-1 models extend the scalar sector of the SM into three SU(3)l scalar triplets: 
one heavy triplet field with vacuum expectation value (VEV) at high energy scale (x) = v x -> 
which produces the breaking of the symmetry SU{3)l®U{1)x to the SM electroweak group 
SU(2)l <8) U(1)y, and two lighter triplets with VEVs at the electroweak scale (p) = v p and 
(77) = v v , which induce the electroweak breakdown. In the version without exotic charges, 
after the spontaneous breaking of the gauge symmetry and rotations into mass eigenstates, 
the model contains 4 massive charged Higgs (H^ 1 , H^), one neutral CP-odd Higgs (^4°), 
3 neutral CP-even Higgs (h°, H®, H®)-, and one complex neutral Higgs {H®) bosons. In 
particular, the charged sector is composed of two types of Higgs bosons: hyper charge-one 
Higgs bosons Hf which exhibit tree-level couplings with the SM particles, and hypercharge- 
two Higgs Hf bosons which show couplings with the SM matter through mixing with non-SM 
particles. An study of production of the hypercharge-one Hf bosons at hadron colliders was 



Production of multiple charged Higgs bosons in 3-3-1 models 



3 



performed by authors in ref. [21] in the framework of the 3-3-1 models. 

In this work we show that the above 3-3-1 model can possess an specialized 2HDM type 
III at low energy, where an special basis is preferred, exhibiting specific couplings in the 
Yukawa constants [22]. Thus, like the 2HDM-III, the 3-3-1 model can predict huge flavor 
changing neutral currents (FCNC) and CP-violating effects, which are severely suppressed 
by experimental data at electroweak scales. One way to remove these effects, is by imposing 
discrete symmetries, obtaining two types of 3-3-1 models (type-I and -II models), which at 
low energy exhibit the same Yukawa interactions as the 2HDM type I and II. In the first 
case, one Higgs electroweak triplet (for example, p) provides masses to the phenomenological 
up- and down-type quarks, simultaneously. In the type-II model, one Higgs triplet (p) gives 
masses to the up-type quarks and the other triplet (jf) to the down-type quarks. One way 
to distinguish a 2HDM from an underlying structure (3-3-1 triplet structures in our case) 
is by identifying multiple like-charged Higgs resonances, for example the two Hf 2 charged 
bosons. To explore this possibility, we show by the method of recursive expansion [23J that 
if mixing couplings with the heavy quark sector of the 3-3-1 model is considered, then the 
hipercharge-two Higgs bosons U 2 can couple with the light SM quarks at tree-level. Thus, 
the H 2 Higgs bosons can be produced in the same production modes as the Hf Higgs, but at 
different statistical rates. In particular, we study the if^-boson production at CERN-LHC in 
pp H2H2 pair production and associated tbH^ single production in the framework of the 
type-I and -II 3-3-1 models. For comparison purposes, we include the if^-boson production 
for the above modes. We also discuss H 12 decay processes to tau leptons (r) assuming that 
— > tv) = 100%. We compare transverse mass distributions for different mass values 
of the Higgs bosons. 

This paper is organized as follows. Section II is devoted to summarize the spectrum 
and interactions of the 2HDM and 3-3-1 models, in particular the interactions with charged 
Higgs bosons. We also show that, at low energy, the spectrum of the 3-3-1 model indeed 
can be separated into a light scale corresponding to the 2HDM type III, and a heavy scale 
decoupled from the light sector. In Section III we perfom the rotations into mass eigenstates 
of the quark sector taking into account small mixing terms between light and heavy scales. 
In Section IV, pair and associated production for charged Higgs bosons are studied in pp 
collisions at LHC collider. Decay into tau leptons is considered in Sec. V. Finally in Sec. 
VI, we summarize our conclusions. 



2 The model 

In order to identify a 2HDM structure into the 3-3-1 model, we first show a brief summary 
of models with two Higgs doublets. 



2.1 2HDM-III Model 

The Yukawa's Lagrangian for the 2HDM type III is as follows 
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+ F L ($ 1 hl ij + ® 2 ht ij )e j R + h.c, (1) 

where q l L = (t/*,/}*)^ are the up- and down-type quark doublets for each left-handed flavor 
component U l : (u,c,t) and D l : (d,s,b) for i = 1,2,3, while 1 % L = (V,e l )i are the neutral 
and charged leptons for three flavors. The coupling constants h ni A are the components of 
non-diagonal 3x3 matrices in the flavor space for f = U, D,e and n = 1,2. $12 are identical 
hypercharge-one Higgs doublets and $1,2 = icr 2 $i 2 are conjugate fields with <j 2 the Pauli 
matrix. Assuming that both doublets acquire vacuum expectation values (VEV), ($1)0 = v\ 
and ($2)0 = V2, the doublets are written as follows: 



The mass eigenstates are related to weak eigenstates by: 




R,,« = ( °t S J' a ) , (4) 

where tan/3 = V2/V1, and a is a mixing angle expressed as a function of parameters from the 
Higgs potential. In the most general 2HDM, the parameter tan/3 is considered as unphysical, 
thus, it is necessary to redefine the Higgs doublets into a basis independent from /3 |24J. 
However, if some type of restriction (discrete symmetries, symmetries from bigger gauge 
models, etc,) is imposed, this parameter can be defined in an special basis, thus it will be a 
physical parameter. 

After the spontaneous symmetry breaking, the Yukawa Lagrangian in Eq. gives mass 
matrices in the quark sector of the form 



M v = —(hV Vl + h 2 J v 2 ), 

M D = -L(h? Vl + h°v 2 ). (5) 

The quark mass eigenstates can be obtained by unitary transformations of the left- and 
right-handed weak eigenstates: U' LR = V^gU^R, D' L R = V® R D L R , where the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix is defined by k = V^V^ . Thus, the matrices in Eq. 
([3]) are diagonalized by the following bi-unitary transformations: 
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mu = ±V^(h^v l + h^v 2 )VE, 
m D = -L V ^{h?v 1 + h°v 2 )V°, 



(6) 



where the rotated Yukawa couplings hf 2 = V^h\ 2 V^ and h'f 2 = V^h? 2 Vg are not in 
general simultaneously diagonalized, leading to FCNC at tree level. However, it is possible to 
suppress these FCNC terms by demanding discrete symmetries in the Lagrangians, obtaining 
two different models: 

Type I: Here, the couplings with one Higgs doublet (for example $ 2 ) is removed from 
the Lagrangian, then only one Higgs doublet (<&i) provides masses to the up- and down- type 
quarks, simultaneously. 

Type IT. in this case the Yukawa couplings are restricted so that one Higgs doublet (for 
example, $i) gives masses to the up-type quarks and the other doublet (<E> 2 ) to the down-type 
quarks. 

Specifically, after the rotations to mass eigenstates, the model exhibits the following 
type-I and -II Lagrangians for the charged Higgs Bosons: 



C 2HDM-I = JJT 



-\/2m D 



K , 



r 2HDM-II 
1~ H ± 



U' 



V 



v 



cot/3 



D' R H + + D' L 



D' R H + + D' L 




where v = \fv\ + v\ = 246 GeV is the electroweak VEV. 



2.2 3-3-1 Model 

2.2.1 Physical spectrum 

We consider a 3-3-1 model where the electric charge is defined by: 



Q = T 3 -^=T 8 + X, (8) 

with T 3 = ^Diag{\, —1, 0) and T 8 = (-^)Diag(l, 1, —2). In order to avoid chiral anomalies, 
the model introduces in the fermionic sector the following (SU(3) C , SU(3)l, U(1)x) left- and 
right-handed representations: 
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Ui : (3*, 1,2/3) 
Q\ = : (3,3,1/3), <! D l R : (3*, 1, —1/3) 

Tr '■ (3*, 1, 2/3) 
Dj 3 : (3*, 1,-1/3) 

l2,3 _ | rr2,3 1 . /Q Q* n\ J 7T 2 ,3 




Qr = \U 23 \ :(3,3*,0),<( (3*, 1,2/3) 

J 2 / : (3*, 1, -1/3) 

l - f ,.1.2.3 1 ./ -| •:> 1 /QA / e i? 2 ' 3 : (!j 1) _1 ) 



where £7^ and .D^ for z = 1,2,3 are three up- and down- type quark components in the flavor 
basis, while v\ and e l L are the neutral and charged lepton families. The right-handed sector 
transforms as singlets under SU(3)l with U(l)x quantum numbers equal to the electric 
charges. In addition, we see that the model introduces heavy fermions with the following 
properties: a single flavor quark T 1 with electric charge 2/3, two flavor quarks J 2,3 with 
charge —1/3, three neutral Majorana leptons {i /1,2 ' 3 )1 and three right-handed Majorana 

12 3 

leptons N R ' ' . On the other hand, the scalar sector introduces one triplet field with VEV 
(x)o = v x , which provides the masses to the new heavy fermions, and two triplets with VEVs 
(p)o = v p and (77)0 = v v , which give masses to the SM fermions at the electroweak scale. 
The (SU(3)l, U(1)x) group structure of the scalar fields are: 



Xi 

\ \ X2 \ (3.-1/3) 

p\ 

75K + 0±<p) ] : (3,2/3) 

TiK + ^^'C?) 

ife I : ( 3 ' -V3). (10) 

After the symmetry breaking, it is found that the mass eigenstates are related to the 
weak states in the scalar sector by [HJ fTTj : 



with: 
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R(3T ' aT ~ C%) ' R ~ ( l) ' (13) 

where tan/3 T = v v /v p , and is a mixing angle obtained from the Higgs potential. 

For the boson-vector spectrum, we are just interested in the physical neutral sector that 
corresponds to the photon A, the neutral weak boson Z and a new neutral boson Z', which 
are written in terms of the electroweak SU(3)l ® U(l)x gauge fields as [IB] . [17] : 



^ = -^i-l(T w yw^ + ^=T w B„ 

where the Weinberg angle is defined as 

S w = Sin(9 w ) = /f 9x - , (14) 
with g L and gx the coupling constants of the groups SU(3)l and U(l)x, respectively. 



2.2.2 Couplings 

Using the fermionic content from Eq. OH]), the neutral gauge interactions for SM quarks reads 
[IT]: 

C N D C = eQMl + t£~9 (9l - + 1 , ffl - g» 76 ) Z'"] g, (15) 

where q isU = (U l , U 2 , U 3 ) or D = (D 1 , D 2 , D 3 ) for up- and down-type quarks, respectively, 
and Q q the electric charge in units of the positron charge e. The vector and axial-vector 
couplings of the Z and Z' bosons are written in table [I] for each SM quark [25J . 

On the other hand, from the kinetic term of the Higgs Lagrangian, we obtain the following 
Higgs-Higgs- Vector interaction associated with the charged Higgs sector: 



iC HHV = - ie [ H +H^ + HtHz]{ P -qYA, 



— [C 2W H+H^ + 2S 2 W H+H^] (p - qf Z, 
+ l(Cwr + 7%) H+H^ + 2(1 + I* ) H+H 2 ] (p - qf (16) 
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Fermion 


9l 


9l 


~9l 


9l 


D 1 


_I _L 2 o2 
2 ^ 3 J W 


l 

2 




-l 




2X/3-4S2, 


D m 


_I _l 2 o2 
2 + 3 J W 


1 
2 








2x/3-4S^ 


u 1 


1 4 Q2 

2 3°W 


1 
2 


— 3— 25^, 


-1+2S2, 




2x/3-4S^ 


u m 


1 4 q2 

2 3°W 


1 
2 


3-85^, 


1 




2x/3-45^ 



Table 1: Vector and Axial couplings of SM quarks and Neutral Gauge Bosons. The index 
m = 2, 3 labels the 3* multiplets. 

Finally, we obtain the following SU(3)l <8> U(l)x renormalizable Yukawa Lagrangian for 
the quark sector: 

+ QiP^i m ^ + QlWii + X^n) T A 

+ T L {v*h J vnm + X *h J xnm )J^ + T L P*h T pnl T 1 R + h.c, (17) 

where n = 2, 3 is the index that labels the second and third quark triplet shown in Eq. ([9]), 
and hfaj are the i, j components of non-diagonal matrices in the flavor space associated with 
each scalar triplet : 77, p, \- 

2.2.3 Yukawa couplings at low energy 

We require the breakdown SU (3)iX(/ {l) x —> SU (2) l X U(l)y in the flavor sector. In order 
to identify the particle content at low energy, we use the branching rules shown in Tab. |2} 
where we identify the following (SU(2)l, U{1)y) left-handed (SM) doublet representations: 

?I' 2 ' 3 = (^w) : (2,1/3), 

/,,1,2,3\ 

If'" = (y A sJ :(2,-l), (18) 

while Tl, J^f and (^ 1,2 ' 3 )i are SU(2)l singlets (which we will call non-SM fermions). The 
right-handed fermions are decomposed into SU(2)l singlets with weak hypercharge equal to 
the electric charge. The scalar sector contains the following hyper charge-one subdoublets: 
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(SU(3) L ,U(l)x) -> (SU(2) L ,U(1) Y ) 



Qi:(3,|)- 


Mi: (2,| 




:(3*,0)- 


v ■ 2,3 
— ► l^2?L 


:(2,I) + Jf 


4' 2 ' 3 = (3, f ) 


,1,2,3 


(2, -1) + (^.3)c . (1> o) 


X:(3,^) — 


>$ 3 :(2,- 


■l) + xo:(l,0) 


i?:(3,^)— ) 


► $2 : (2, - 


1) + 4 : (1, 0) 



p:(3,|)-^$ 1 :(2,l)+p+:(l,2) 



Table 2: Branching rules for the (SU(3)l, U{\)x) — > (SU(2)l,U(1)y) symmetry breaking. 



1 " (^K + 1 + <p)) :(2 ' 1} 

tiK+^-^a . (2 _ 1} 

V2 J 

*3=Q):(2,1), (19) 

while xo = (Vv / 2)('Cx ^ ^Cx) an d ^3 are hip er charge-zero singlets, and pf (which according 
with (1121) are identified with the charged Higgs bosons) are hipercharge-two singlets. In 
(fl9"|) we define the conjugate field as $ 2 = icr^%- I n ^ ne above decompositions, the U(l)y 
weak hypercharge, the U(l)x charge and the electric charge were related by: 

Q = T 3 + Y/2 = T 3 - -±=T 8 + X. (20) 

V3 

Taking into account the above branching rules, the Yukawa Lagrangian in Eq. ( 1T7|) can 
be separated as follows: 




+ 7 L {^Kj + ^Kj + ^ h xv) D R 

- C Y (T\j n ), (21) 

where £y(T x , J n ) are mixing terms with the non-SM fermions. The Lagrangian in (I2T!) 
exhibits the same form as the 2HDM-III in Eq. (CQ) for the quark sector, plus extra mixing 
terms associated with the heavy scalar doublet $3 and the other non-SM particles. Thus, for 
v x 3> v PjV , the 3-3-1 model exhibits an effective 2HDM-III couplings in the quark sector at 
low energy. Furthermore, due to the nonuniversal U(l)x values exhibited by the spectrum 
in (Q and fflOl) . not all couplings between quarks and scalars are allowed by the gauge 
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symmetry, which leads us to zero-texture of the Yukawa constants in Eq. ( 1211) |22j . If we 
consider a low energy scenario in which the particles at the heavy scale decouple from those 
at light scales, the quark mass eigenstates may be obtained separately for each scale by 
unitary transformations of the left- and right-handed weak eigenstates: U' L R = Vl R Ul,r, 
D' L R = V^Dl^r, and J' L R = V^ r Jl,r, while the single flavor T 1 quark decouple from other 
components, obtaining that T'£ R = T\ R . In this limit, after the symmetry breaking, and 
using the VEVs from Eq. (fTUj) . we obtain the following mass terms for the SM quarks: 

M v = -^=(h u p v p + h^v v ), 

M D = -L(h°v p + h°v v ), (22) 

which exhibit the same form as the matrices in Eq. (JH]) but with the change V\p — > v p ^ and 
ft<i,2 — ► hp,ri- Then, using bi-unitary transformation analogous to (j6]), we obtain the diagonal 
mass matrices 



m D = V^MnVg = {h D p v p + h°v v ) Vg. (23) 

By requiring appropriate discrete symmetries, we can restrict the Higgs couplings in the 
Yukawa Lagrangian. In particular, we may generate type-I and -II models analogous to the 
2HDM-I or -II if we require the following generalized discrete symmetries on the scalar and 
right-handed fermion fields: 



Type I: 



p 


— >■ 


p, v->-v> x -> 


-x 


Ur 


-)■ 


Ur, Dr -> Dr 




Tr 


— > 


—Tr, Jr — > —Jr 




Nr 


— >■ 


-N r , e R ->■ e R 




P 


->■ 


p, T] ->• -77, X ->■ 


-x 


Ur 


->■ 


E/ fl , -Dr -»- -D R 




Tr 


—> 


—Tr, Jr — >■ — J/? 




Nr 


—> 


-iVfl, ->■ e fl . 





Type II: 



(24) 

In type-I models, one Higgs triplet (p) provides masses to both the up- and down-type 
quarks, while in type-II models the triplets p and rj give masses to the up- and down-type 
quarks, respectively. Thus, eqs. (1231 becomes: 



Type I: m um = -±V?™ 

Type II: m „ m = -^Vf <"» (<<>*)) V^'. (25) 
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If z/ P)7 j = Ux t 2, the mass eigenstates in ([3]) are the same as the mass eigenstates in Eq. ( TTTj) . 
In particular, we can identify the hypercharge-one Higgs bosons Hf as 2HDM-like charged 
bosons, while the other charged Higgs bosons (hipercharge-two) Hf are new charged scalar 
particles beyond the 2HDM. 



3 Yukawa Lagrangian with mixing couplings 

If we consider the complete Lagrangian in (|T7|) (including couplings with non-SM fermions), 
we find the following mass terms after the symmetry breaking [22] : 



(C Y ) = faM^R + h.c. = (U L ,T L ) M UT (^A + (pl,~h) M DJ + h.c, 



(26) 



where Ul,r = (U 1 , U 2 , U 3 ) l,r are the left- and right-handed up-type quark flavor vectors, 
Dl,r = (D l , D 2 , D 3 ) l,r the corresponding down-type quark vectors, Jl,r = (J 2 ,J 3 )l,r are 
two-dimensional vectors associated with the heavy quarks with electric charge —1/3 in (J9j) 
and Ti t R is the single component of the heavy quark with charge 2/3. The mass matrices 
have the following structures in the basis (U,T) and (D, J), respectively: 



M, 



Mught flight \ . ^ 

Gheavy ^-heavy J 



Mi 



UT 



M, 



D.I 





k 


K 


M T 


M D 


s 


S 


Mj 



(27) 



where Mjj, k, K, and are 3 x 3, 3 x 1, 1 x 3, and lxl matrices, respectively, while Md, 
s, S, and Mj are 3 x 3, 3 x 2, 2 x 3, and 2x2 matrices, respectively. The components of 
the above mass matrices are: 



Mlight 

^-heavy 
flight 
Gheavy 



T(J) 



V2 



k(s) 



V2 



K{S) = ^h^v x , 



(28) 



where M Hght and fu ght are of the order of v p , v v ~ GeV, while A heavy and G heavy are of the 
order v x 3> GeV. The diagonalization of the matrices in ( 1271) can be obtained by unitary 
transformations of the left- and right-handed weak eigenstates: 



tP'l,r = 0\r^l,r, 



(29) 



12 



R. Martinez and F. Ochoa 



where ip' L R = (U', T')l,r or (£>', J')l,r are the mass eigenstates. In order to find the form of 
the matrices Ol,r, we separate them into two rotations: 



O 



L.R — Ml,rWl,R — ( 




-B 



L,R 



Bl,r 

1 



(30) 



that diagonalize the 
V[ R or 1 that rotate the J or T 



where Vl,r = V^'r are ^ ne same bi-unitary transformations as in 
U and D components in the low energy limit, and Pl,r 
components, respectively Since the first rotation through Ul,r does not lead to a completely 
diagonal mass matrix due to the mixing terms fu g ht and Gheavy in (|27|) . we choose bi-unitary 
rotations B^ R = B L ' R by requiring the vanishing of the off-diagonal components of the 
following matrix: 



0\M^O R = WMmMrWr 



( might o 
V U ITlfieavy 



(31) 



Considering the scenario of small mixing terms (near the decoupling limit) (Gheavy 
A-heavy, flight "C Mught) and using the method of recursive expansion [23J, the authors in ref. 
shows the following solutions: 



B L « (m,Gt + M) (m 2 h ) 1 
r ~ G(m ft ) _1 , 



5 



(32) 



where m ; = V^Mh^Vr, m h = PlA heavy P R , f = Vlfi ight PR and G = P\G heavy V R . Replacing 
fl32|) in fl30|) and (129]) . we obtain the following mass eigenstates: 



R 



T'r 



(V^)U L - {B U L )T L » (VnU L 



(VEB U L yU L + T L : 
(V^)Ur - (B R )T R « (T^C/* 



/ / brp 



t4 + T L 



k 

in? 
U R + T R 



-R 



(33) 
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D' 



D', 



R 



{V?)D L -{BZVp)J L *{V?)D L 



m 



J'l = (V L D B^D L + (V^)J L 



m D 



m 



mj 



Di 



Jl 



—V J] 
mj R 



J, 



R 



J'r = (VgB°y DR + (v£)J R 



—V D] 
mj R 



Dr + (V£)J r . 



(34) 



With the rotations from (I3"3"|) and (I3~4j) for the quarks, and (TTT|) - (ir2|) for the Higgs, the 
Yukawa Lagrangian of the 3-3-1 model in (JT7J) can be written in mass eigenstates. In partic- 
ular, for the couplings between the charged Higgs bosons and the (U', D') mass eigenstates, 
we obtain: 



- & = mv?) [(-hfSft, + h°C PT ) Vg - (~h J p Sp T + h J v C, T ) V£B$ 
+ WV^) [(-h u p S, T + h"Cp T ) V" - (-h T p S, T + h T v C, T ) B U R [ 

' r h D p Vg + h}B$ 



D' R Ht 



+ mBM ] ) [~h U p VK + h T p B u R [ 
+ U[{B U L ') \-h D p Vg + h J p B u R '} D' R Hi + h. 



.c. 



(35) 



Since the rotations B^' R are suppressed by the inverse of the masses of the heavy T and 
J quarks according to ( l32i) (where nih = Tn^j), the couplings associated with h J '^ in ( 1351) 
are negligible with respect to the h p ^ couplings. In addition, other terms can be removed 
by the discrete symmetries in (1241 . Taking into account the above facts and using Eqs. (|25|) . 
the Lagrangian (1351) for type-I and -II models can be written as: 



L 'Y 



= U' L 

+ d£ 

+ h.c. 



■\/2m D 



K 



V 

■\r2mu 



T, 



Pt 



D' R Ht + D' L 



T, 



v 



Pt 



vCp 7 



U' R H^ + U' T . 



^ mD {BVvE) 



vC. 



Pt 



D' R H+ 



(36) 



MI 



+ D> 



y/2mr>K 



v 



cot fa 



D' R Ht + D' L 



T, 



v 



Pt 



U' R H^ 



^(B?V T J W, 



vC. 



Pt 



L " L L 



U' R H^ + h.c. 



(37) 
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where k is the CKM matrix and v = \/v 2 p + v 2 , = 246 GeV is the electroweak VEV. We see 
that the couplings of the Hf bosons in the above Lagrangians are analogous to the if ± -boson 
couplings of the 2HDM in (jTJ). Thus, we identify the Hf bosons with a 2HDM-like charged 
Higgs bosons (furthermore, if at low energy we neglect the coupling with the Z' boson in 
Eq. (USD, we obtain that Hf ->■ 



4 Production of charged Higgs bosons 

Since the couplings found in (|36|) and (1371) are proportional to the quark masses, the largest 
contribution comes from the top quark (m t ~ 173 GeV). Thus, considering the couplings to 
the third up type family, the dominant contribution of the Yukawa interactions in (136!) and 
flHID are: 



where we consider K t b ~ 1 for the (3, 3) component of the CKM matrix. We are interested 
in comparing the production ratios between Hf and Hf, then in the above Lagrangian we 
also assume for simplicity that V^'^ ,J ~ 1 for the U, D and J quarks. It is to note that 
type-I and type-II models exhibit different Yukawa couplings in the down sector, but the 
same couplings in the up sector. Since we neglected the contribution from the down-type 
masses, we see in Eq. (|38|) that both models possess identical couplings. 

Taking into account the couplings in ( fl5|) . ( fTBl) and ( 138]) for type-I and -II models, we 
may study the production of Hf 2 Higgs bosons. Fig. [1] shows different partonic processes for 
Higgs production in pp collisions, where (a) corresponds to quark-antiquark annihilation for 
pair production through vector electroweak bosons, while (b)-(f) corresponds to associated 
production tbH 12 through gluon-gluon collisons ((b)-(d)) and quark-antiquark annihilation 
((e) and (/)). 



4.1 Pair production 

For the pair production from Fig. QJa), we use the couplings in (TIB"]) . We observe that 
Tp T , M Hl 2 and M z > are free parameters, while and gx can be parametrized as func- 
tions of the electric charge of the proton e and the Weinberg angle: e = giSw and 
gx = V^giSw/ '\/3 — 4:S%y (from definition in (T141). The parameter Tp T only appears in 
the coupling of Hf with Z' through the coefficient {C2p T + Tyy) in (TIB"]) . Just for numerical 
purposes we take Tp T = 10. At CM energy of 14 TeV in pp collisions, we use the Calchep 
package |26j to obtain the pair production cross sections for Hf and Hf bosons. Fig. |2^a) 
shows the dependence of the cross section with the mass of the Z ! boson for three different 
Higgs masses: M Hl 2 = 300, 500 and 1000 GeV. The behaviour shown by the curves can be 
understood as follows: For M z > < 2M Hl 2 , the production of real Z' bosons is supressed by 
kinematical conditions, so that the contributions come only from A and Z bosons. Above 
the resonance region Mz> = 2M# 12 , the production of Z' increases, which leads to larger 
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cross sections. However, when the Z' boson becomes heavier, the energy of the collision 
is not large enough to obtain an appreciable Z' production, thus the A and Z contribu- 
tions are again the dominant modes for pair production. On the other hand, we see that 
around the region of the Z'-resonance, the cross sections splits in two branches, where the 
H2 bosons exhibit larger contributions than the 2HDM-like bosons Hf. This splits is due 
to the different contributions exhibits by the Hf and bosons with the gauge bosons in 
Eqn. ffl6|) . Using the notation gUiHiV to designate the Higgs-Higgs- Vector couplings, and 
with SM numerical inputs, we obtain the following relations: 

9h 2 h 2 a 
Qh^a 
9h 2 h 2 z 
Qh^z 
9h 2 h 2 z' 

g^HrZ' 

which indicates that due to the Z'-couplings, the ratios between cross sections a(H^)/a(H^) 
can be as large as 3.7 2 ~ 13.7. i.e as one order of magnitude larger. Thus, we can see in Fig. 

for the curves M Hl2 = 1000 GeV, that the boson exhibits a peak which is about 
one order of magnitude larger than the Hf peak. 

On the other hand, Fig. Wib) compares the dependence of the cross sections with the 
Higgs masses for Hf, H 2 1 an d H 2HDM Higgs bosons. We use the value Mz> = 2 TeV. As 
with the Fig. EJja), the cross sections splits into various branches for Mh 12 < Mz>/2 = 1 
TeV, where the largest contribution comes from the H 2 bosons. In particular, we see that 
for 2HDM we obtain smaller production ratios due to the fact that in this model the Z' 
contribution does not exist. 

4.2 Single production 

Figs. W[b)-(f) corresponds to single Higgs production in association with bt quarks (we 
show the production of the iff 2 component. The Hj~ 2 is identical but in association with bt 
quarks). In this case, in addition to (fl6l) . we must take into account the parameters from 
couplings in ( fl5l) and ( )38l) . Thus, we have the new free parameter B® for the quark couplings 
with H 2 . This parameter corresponds to the mixing rotations associated with the down- type 
mass matrices Mpj from fT27j) . According to Eq. fl32l) . and for the third Down-type family, 
we obtain the following relation: 

Bl w (m b S j + sm^ (m^y 1 (40) 

where mj, rj 4 GeV corresponds to the mass of the b quark, mj ~ v x 3> GeV is the mass 
of the heavy J quarks, S = S ~ v x the heavy mixing component of Mdj in (1271) . and 
s = s ~ v p ^ v ~ GeV the corresponding light mixing components. Due to the above order 
of magnitudes of the VEVs, the mixing rotations in fT40|) exhibits small values (B® C 1). 



( ~'2W 

(C 2 8 T + T w ) 
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However, the H 2 couplings in Eq. ( 1551) depend on the ratio B£ /Cp T . Thus, although the 
mixing rotations takes small values, the couplings with the Higgs H 2 may be appreciable 
for small Cp T . Then, we use the parameter r = B® /Cp T instead of Bj? for our analysis. 
First, we obtain the cross section for associated production of Hf Higgs bosons. Fig. |5Ja) 
exhibits the dependence of the cross section with the parameter Tp T for the three Hf masses, 
where large Tp T values lead to large production cross sections. Second, for the if^-boson 
production, we calculate the dependence of the cross section with the parameter r, as shown 
in Fig. [5J6). Since the Yukawa couplings is proportional to r, the cross sections increase 
with r. 

4.3 H^ 2 tb production 

In the above subsection we obtain the cross sections for single production in association 
with bt quarks. However, we can obtain the same final states from pair production in the 
diagram [T^a) if we suppose that one of the Higgs bosons decay to bt. Thus, the total cross 
sections a(pp — > tbHf 2 ) have contributions from all the modes shown in Fig. [TJ To explore 
how large is the contribution from each mode, we calculate the corresponding cross sections 
using the following values: M z > = 2 TeV, Tp T = 10 and r = 0.1 (we note that these values 
are consistent with the condition of Bj? small: Bj? = rCp T ~ 0.01 <C 1). Figs. 01(a) and 
(6) display the dependence of the cross sections with the Higgs masses for Hf and H±, 
respectively. We see that, in general, the dominant modes come from Gluos-Gluon collisions 
(Figs. \V[b)-(d)). We also see that, as it is to expect, the contribution from Higgs pair 
production qq — )■ Z' — >■ H^ 2 H^ 2 tbHi 2 , becomes resonant for Mh x 2 = Mz>/2 = 500 GeV. 
Furthermore, this mode becomes larger than the gluon collision for Hf masses around the 
resonant region (see Fig. H3^)). The mode from Fig. Wif) is, on the contrary, negligible for 
all the range of Mh 12 - 

5 Decay of the charged Higgs bosons to leptons 

To explore possibilities to distinguish different like-charged Higgs bosons, we study the decay 
mode Hf 2 — > ru from associated Hf 2 bt production. To have this decay, we first must study 
the lepton couplings in the Yukawa Lagrangian. From the spectrum in (Q and (1T0|) . assuming 
Majorana mass terms for the singlets Nr, and according with the discrete symmetries in (I24p . 
we obtain the following renormalizable Yukawa Lagrangian for the lepton sector: 

+ \ E ?(if) c V(/ a ^^ fi ,«+k, (4i) 

a,/3,7=l 

where i, j = 1, 2, 3 label the family index (three Ll triplets), a, /3, j — 1, 2, 3 label the flavor 
components into each triplet, and e°^ 7 corresponds to the antisymmetric Levi-Civita tensor 
(1 for even permutations of (123), —1 for odd permutations, and if any index is repeated). 
For the charged sector, we obtain the following mass matrix after the symmetry breaking: 
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M e = —h;v p . (42) 

For the neutral leptons, the Lagrangian in (j4"T|) contains mixing terms which produce neutrino 
mass terms. The authors in ref. [27] study three see-saw scenarios in order to obtain small 
neutrino masses. On the other hand, as with the quark Lagrangian, we rotate the Higgs 
sector to mass eigenstates according to Eqs. (ITT)) and f lT2|) . In particular, for the couplings 
with the charged Higgs bosons and the SM leptons (e\ z/), we obtain: 



- £?' 2 = vl (-h^Sfr) e^Ht + el (-h p ) vgHj + h.c. (43) 

Using Eq. ( 142]) . we can write the coupling constant h e p in terms of the lepton masses, 
while v p = vC/3 T . Thus, for the tau lepton (r), Eq. (143l reads: 



_ £ - M = V Pt — TrH + _ hpTEu C RH - + h (44) 

In order to obtain small neutrino masses, the see-saw mechanism in 3-3-1 models requires 
small values of the parameter h p (typically h p ~ 10 -7 — 10~ 4 ). However, in the framework of 
the inverse see-saw mechanism, it is possible to obtain couplings as large as h p ~ lGeV/v p 
which is consistent with small neutrino masses [27) . 

With the above considerations, we calculate the cross section distribution in pp colisions 
in the CERN LHC hadron collider, based on an integrated luminosity L = 100 fb _1 at CM 
energy of 14 TeV. We assume that Br(Hf 2 —t tv) = 100%. We also consider the dominant 
contribution from gluon-gluon collisions (Figs. [[](&)- (<f)), and keep large Z' masses to neglect 
its effects. We fix the following parameters: 



Tp T = 10 
r = 0.1 

hp ~ lGeV/v p = \GeV/vCp T w 0.041. (45) 

Fig. [5] shows the Higgs boson tranverse mass distributions for the r final state. The 
distribution in Fig. [5](a) shows the Higgs events for Mh-l = 300 GeV, and for Mh 2 = 500 
and 1200 GeV. For bosons with mass 500 GeV, the signal is overlapped by the Hf 
background, while for 1200 GeV bosons we see a small excess over the background. Fig. 
\Mb) shows the same distributions but for Mh x = 500 GeV, which overlaps the signals. 
Finally, Fig. [5](c) displays distributions for Mjj 1 = 1000 GeV, where an observable 
signal with mass 500 GeV arises over the H l background. Although the production of 
the hipercharge-two Higgs boson is in general small in relation with the hyper charge-one 
Higgs, we see from the above analysis possible scenarios where the two signals may be 
distinguishable. Furthermore, this small signal could be improved if more sophisticated 
discriminating distributions are used, or other decay channels are considered [28j . 
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6 Conclusions 

The identification of multiple Higgs boson signals could reveal many features about the 
underlying model beyond the SM. For example, if a charged Higgs boson is detected, further 
analysis will be necessary to test the compatibility of different models with the experimental 
data. In this paper we show that both versions of the 2HDM (type-I and type-II models), 
which contains one charged Higgs boson, may be embedded into a 3-3-1 model, which exhibits 
two charged Higgs bosons: a hipercharge-one Hf and a hyper charge-two H 2 Higgs boson. 
At low energy, Hf can be identified with the charged Higgs boson of the 2HDM, while H 2 
are other bosons from the underlying 3-3-1 model. Thus, in this case, the identification of 
two like-charged Higgs boson signals may reveal new physics beyond 2HDM. Taking into 
account mixing couplings between the two scales of the model [v ~GeV and v x 3>GeV), we 
show that the H 2 can be produced through the same production channels as the H 1 . Using 
the method of recursive expansion, we found that after rotations to mass eigenstates, Yukawa 
couplings between the SM fermions and H 2 arise due to the small mixing angle associated 
to the rotation matrix B£. Furthermore, since the Yukawa couplings appear through the 
ratio /Cp T for type-I and -II models, events of H 2 may be enhanced to observable scales 
if Cp T takes small values. We show that pair production of charged Higgs bosons can be 
significantly enhanced due to the contribution of a heavy Z' neutral gauge boson, predicted 
by the model. The dominant mode for associated btHf 2 is through gluon-gluon collisions. 
However, pair production through quark anhilation can be as large as single production due 
to resonant intercharge of the heavy Z' boson in 3-3-1 models. By considering decays to 
leptons H 12 — > rv T1 we obtain scenarios where small peaks of -boson events in transverse 
mass distributions can be identified over the Hf background. This small signal could be 
improved if more sophisticated discriminating distributions are used. Thus, in case that 
charged Higgs bosons are detected, the identification of multiple like-charged Higgs boson 
signals may be a possible discriminating method to test different theoretical models beyond 
the SM. 
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Figure 1: Charged Higgs boson production in pp collisions for (a) pair production in quark- 
antiquark annihilation, (b)-(d) associated single production in gluon-gluon collision and (e) — 
(/) quark-antiquark annihilation 
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Figure 2: Dependence of the cross section for charged Higgs-boson pair production with 
(a) the Z'-boson mass for M H = 300, 500 and 1000 GeV, and with (b) the Higgs mass for 
M z > = 2 TeV in 3-3-1 and 2HDM models. 




Figure 3: Dependence of the cross section for associated (a) tbH l production with tan/3^ 
and (b) tbH^ production with the ratio r = B%/C Pt for M Hl 2 = 300, 500 and 1000 GeV 
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Figure 4: Dependence of the cross section for (a) t&i^ and (6) t6iJ 2 productions with the 
Higgs mass in different production modes. We use the inputs tan/3y = 10, r = 0.1 and 
M z , = 2 TeV. 



Production of multiple charged Higgs bosons in 3-3-1 models 



25 





200 400 600 800 1000 1200 1400 

mj, (GeV) 

(c) 



Figure 5: Transverse mass distributions in H± 2 — > tv t decay in associated Higgs-boson 
production for (a) M Hl = 300 GeV, (b) M Hl = 500 GeV, and (c) M Hl = 1000 GeV compared 
with Hf bosons for M H2 = 500 and 1200 GeV. 



